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1. Introduction 

This short paper can be thought of as an extension of [1] where we have proved 
that a regular real model set with a real internal space is, up to translation and 
alterations of density zero, uniquely determined by its 2- and 3-point correlation 
measures. There we have also given an example in the case that the internal space 
is the product space of a real space and a finite group that shows that this result 
does not extend in general to more complicated internal spaces. More precisely, 
we have given an example of two distinct model sets having equal 2- and 3-point 
correlation measures, created using a single cut and project scheme, that are not 
translationally equivalent, even allowing for alterations of density zero. 

In this paper, extending the setting to weighted model sets, we offer examples of 
pairs of weighted model sets which have equal 2, 3, 4, 5-point correlation measures, 
created using a single cut and project scheme, that are not translationally equiva- 
lent. These pairs are created by imposing a 6-colouring on a previously constructed 
aperiodic model set and then further imposing on them two different weighting 
schemes. The resulting weighted model sets are then different, not even having 
any weights in common, but their correlations measures, up to the fifth one are 
identical. The relevance of this type of result to the theory of long range-aperiodic 
order and crystallography and some comments on its history are discussed in 

The key to this result (and also to the previous examples for unweighted model 
sets) is to use corresponding results, developed in [5j, for one dimensional periodic 
sets. There is a simple way in which to intertwine this periodicity with the structure 
of an ordinary model set, and this leads to aperiodic weighted model sets with 
points of several types or colours which have the same coincidence of 2, 3, 4, 5-point 
correlations. 

In brief outline, we start with a weighted periodic model set on the real line 
whose internal space is Z/iVZ for some > 1. We let S — (R'' x H,C) be an 
arbitrary cut and project scheme, where _ff is a locally compact Abelian group and 
£ is a lattice in M*^ x H. As usual, we denote by L the projection of £ into R'^ and 
denote by (•)* the star map of 5 from L to H. Now assume that there is a surjective 
homomorphism a : L — > Z/A^Z. Then there is a combined cut and project scheme 
specified by = (R'* x {H x Z/NZ), £'=),wherc C := {{x, {x\a{x))) -.xeL}. We 
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call a{x) the colour or type oi x L. We use this new cut and project scheme to 
create coloured model sets with several windows, each corresponding to a different 
colour. In order to get a weighted model set we weight the points of this coloured 
model set according to their type. 

We assume that the reader is familiar with the basic theory of model sets [6]. 
In §2 we briefly recall the definitions that we need and the important concepts 
of uniform distribution and the point-correlation measures. Closely connected to 
correlation measures are the pattern frequencies, which are generally more conve- 
nient for our purposes. In §3 we consider discrete periodic point sets on the real 
line, regarding them as model sets, and present the general formula for their finite- 
point correlation measures. In §4 we elaborate the method of building an extented 
weighted model set from a weighted periodic system and an aperiodic model set, 
and then determine the precise form of the resulting pattern frequencies. In §5 we 
take an example from [S] of two weighted periodic systems (based on Z/6Z) which 
have equal 2, 3, 4, 5-point correlation measures (we offer a short proof of this), and 
then use it to show that our extension construction will produce any number of 
aperiodic weighted real model sets with equal 2, 3, 4, 5-point correlation measures. 
Finally, by way of illustration, we offer an example of this using a 6-colouring of 
the vertices of Franz Gahler's shield tiling. 

2. Model sets and correlations 

We work in M.'^. The usual Lebesgue measure will be denoted by £. The open 
cube of side length R centred at is denoted by Cr. We begin with the cut and 
project scheme S = (M'', H, C) consisting of a compactly generated locally compact 
Abelian group H and a lattice C C K'* x H for which the projection mappings tti 
and 7r2 from K*^ x H onto R'' and H are injective and have dense image respectively: 



X ^ X ^ X* . 

Then L := 7Ti{C) is isomorphic as a group to £ and we have the mapping (•)* : 
L — > H , X ^ X 1-^ X*, with dense image, defined by 7r2 o (7ri|£)^^. 

The statement that £ is a lattice is equivalent to saying that it is a discrete 
subgroup of R'' x H and that the quotient group T := (R'' x H)/C is compact. We 
let 6h be a Haar measure on H, normalized so that the product measure £ (8) Oh on 
R'^ X H gives total measure 1 to a fundamental region of the lattice C. We let 9j 
be the canonical Haar measure on T whose total measure is 1. 

For W CH, 

A{W) := {ue L:u* eW}. 

A set W C H is called a window if S° C C S for some compact set S C i? 
which satisfies E° = S. 

We shall assume throughout that all windows that we use have boundaries of 
9 H -measure 0. 



^Othcr normalizations arc possible. The one we have chosen leads to the formula of Thm[T] 
Other normalizations produce multiplicative factors in this formula. None of this is relevant to 
what we need here. 
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We will deal with multiple windows, and it is convenient to allow windows to be 
empty (which is allowed by the definition). 

A (regular) model set or cut and project set is a set of the form A = a; + yl(W) 
where is a window and x G M''. In the sequel we shall only have need of the 
simpler model sets of the form A = yl(M^). 

Model sets are uniformly distributed point sets: 
Theorem 1. [7] Let W C H be a window. Then 

lim -i— card(/l(W^) n Cr) = OniW). 

R—>oo l[L^ R) 

□ 



Let m be any positive integer and let m :— {1, . . . , to}. We assume that we are 
given disjoint windows Wi , . . . , Wm , and then the corresponding model sets 

Aj = A{{Wj) c M''. 

Let A := [J^^ A^ (a disjoint union). We assume given a set w — (wi, . . . , Wm) 
of weights Wj G K. For a; G K'^ we define 



w{x) 



Wj if a; e Aj 
otherwise. 



In our examples below we shall only use non-negative weights. 

We call (Ai, . . .A™) with the weights w a weighted model set. More generally 
one would allow arbitrary translations of the colour component sets A; provided 
the translated point sets do not overlap, but we have no need of this here. For 
notational simplicity we shall use the symbol A'^ (or often simply A is the context 
is clear) to denote the coloured/ weighted model set that we have just described. 

The n+ 1-point correlation {n = 1,2,...) of a model set A'^ (or more generally 
any weighted locally finite subset of M'*) is the measure on (R"*)" defined by 

7a"'^^^(/) = lim TTTTT Yl w{x)w{yi) . . .w{yn)f{-x + yi, . . . , -x A 

^ yu...,yn.,'j:<£CRnA 

= J™ TTTTT Y] w{x)w{yi)...w{yn)f{-x + yi,...,-x + yri), 

^ ' xeCRnA 
yi,...yneA 

for all continuous compactly supported functions / on (R*^)", [3]. The simpler 
second sum is a result of the compactness of the support of if of / and the van 
Hove property of the averaging sequence {Cr} (namely that the measures of the 
if-boundaries of the sets Cr have vanishing relevance relative to the total volume 
of Cr as i? ^ oo . 

Because model sets are Meyer sets, the sets of elements yj — x which make up the 
values of the arguments of / occuring in the sums lie in the uniformly discrete set 
A - A. For any z = {zi, ... ,Zn) G (A - A)" and any fc = (A:(0), . . . , k{n)) G m"+^ 
we can count the occurrences of zi, . . . , z„ in the form —x + yi, . . . , —x + y„ where 
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X e Afe(o) , Vj e Ak(j) , j = l,...,m, namely 
freqfc({zi, . . . ,z„}) = 

i™ err^ x Cardjx € Cr : x € Ak(o) , Vj e Ak(j) ,x = -Zj + yj,j = I, . . . ,m} 

m 

= TTT^^^'^^'f^ e L n Cij : a;* G Wfe(o) n Pi (-4 + W^feo))} 

I VKfe(o)nf|(-z; + M^fe(,)) 

where we have used the imiform distribution theorem. 

Hence for model sets we find that all these correlation measures exist and 



(2) 7^"+^^ = ^ ( ^ ■ ■ ■ fi'eqfe(2;i, ...,Zn) \ (5(^1,., 

zi,...,z„eA-A \feem"+i / 



3. A PERIODIC EXAMPLE 

Let A'' be a positive integer and let H = Z/NZ. Then we have the cut and 
project scheme: 

M R X Z/NZ Z/NZ 

(3) ^ y 

X < > (a;,a;Ar) 1-^ xn , 

where xn ■= x mod N and Z := {{x,xn) ■ x G Z}. Let m be as above. We 
assume that we are given disjoint subsets (windows) Ai, . . . ,Am C Z/NZ and 
corresponding colour weights wi, . . . , Wm- The corresponding model set is A = UAj 
where Aj = A{Aj) = {x € Z : xn G Aj}, j = l,...m. Each colour set A^ is 
periodic, repeating modulo N, containing all those integers congruent mod N to 
an element of Aj. 

For each j = Q, . . . N — 1 define 

(4) Cj 



Wk if j G Ak 
otherwise. 



A given pattern (ri, . . . r„) of integers modulo N must occur in A in the form 

s + ri, . . . , s + r„ mod A, s = 0, 1,...,A — 1. 

If it occurs with colours A:(l), . . . , k{n) and k{0) is the colour of s then rj G A^j) 
Ai.(o) for all j, and the weighting is 

Wfe(o) • ■ • Wk{r) = CsCs+n ■ ■ ■ Cs+r„ • 

Thus, for this weighted model set we have 

7a ^^^ = ;^ X] Mn{ri,...,rn)S^ru-,r,.) 
(ri,...,r„) 
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where the sum runs over all possible patterns of length n and where 

N-l 



s=0 



(5) M„(ri,...,r„) := XI 

By the above arguments, 



(6) Af„(ri, . . . , r„) = ^ w^qj . . . Wfc(„) O^/nz A^q) n f] (-r* + A^^)) 

4. Aperiodic Examples 

Start with the cut and project scheme 5 of [J2l along with the corresponding 
notation. Let 

a: L — > Z/NZ 

be any surjective homoniorphisni with the property that 

{{x*,a{x)) : X e L} is dense in H x Z/NZ . 

Form the new cut and project scheme S'^ — (M'',i7 x Z/NZjC^) where := 
{{x, a{x)) : X e L}: 

R'' < — R"^ X H X Z/NZ — > H X Z/NZ 

a; < — > {x,a{x)) i-^ {x*,a{x)) 

is clearly discrete and since the index \Cf' : {\x G £ : a{x) ~ 0},0)] is finite and 
the group Z/NZ is finite, the quotient of x H x Z/NZ by is compact. In 
short, C is a lattice inM.'^ x H x Z/NZ. 

Now we let m and the sets Aj be as in ^ Let be a non-empty window in H 
and set Wj :— W x Aj d H x Z/NZ. This produces from S"^ coloured model sets 

(8) K'j=A''{Wj). 

Let A*^ := U^'L]^A|. Notice that the actual points of the model sets involved here 
form a subset of the model set determined by the original cut and project scheme 
S, whereas the colours are being determined by the periodic scheme. 
The (n + l)-point correlation for A^ is 

(9) 



2:i,...,z„eA'=-A'= \fcer7x" + i / 



where 



freqfe({zi, . . . , Zn}) = Ohxz/nz I Wk(o) D f] {-z* + Wk(j)) 
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Let z* = {v*,r*) x Z/7VZ. Then for {q, r) e H x Z/NZ, 

n 

{q, r) e Wk(o) n p| {-z* + W^^)) 

^ {q,r)eW X Afe(o) and {q, r) 6 {-v* , -r*) + W x Ak(j) 
<^ qeWn f]{~v* + W) and r G A^o) D f]{-r* + A^^)) . 
Thus the (relative) frequencies are given by 

freqfc(zi, . . . , Zn) = Oh r\ f] (-u* + W)^ 9z/nz ^^fc(o) n f] (-r* + A^^,) 

The first term of this factorization is independent of k and as a consequence we 
can rewrite (ITOl) as 



(10) 7i"^'^ = E Onlwrif] {-V* + W) M„(ri, . . . ,r„) <5(,„...,,„) 

zi,...,2„eA<=-A'= \ i=i / 

It is not particularly important for our purposes that the frequencies here be 
absolute. As we already pointed out, that depends on normalizing the Haar mea- 
sures so that the corresponding Haar measure on (R'' x H x Z/NZ) / has total 
measure equal to 1. What is important is to realize that if we colour and weight by 
a second set of weights , . . . , w'^ and a second set of windows A'^ , . . . , A'^ so that 
the expressions of © are equal for some n, then also the [n + l)-point correlations 
of the two corresponding weighted model sets arising from using one or the 
other of these two colour/weighting schemes will be the same. This information is 
contained in ([5]), and it is this form that we shall see in the examples. 

It is interesting to note here that the formula for the frequencies makes it look as 
if we are dealing with a simple product structure. However, the points Zj = {vj, rj) 
are not truly from an unrestricted product. In fact, rj — a{vj). The reason for the 
frequencies to be given as they are is that the lattice already has this special 
structure built into it. We have assumed that its image is dense in H x Z/NZ, so we 
have a cut and project scheme, and this allows us to use the uniform distribution 
of model sets to derive the frequencies in terms of the measures of the windows. 

5. Examples 
With = 6, there are two sets of weights 

wsi [11,25,42,45,31,14] 

(11) WS2 := [10,21,39,46,35,17] 

which, when used to weight the sets Aj = A'j = {j} mod 6, j = 0, . . . , 5, determine 
identical results on the left side of (O for n = 1, 2, 3, 4, see [5] §5.3. It follows from 
our discussions that the corresponding weighted model sets built in ([51), though 
quite different, nonetheless have equal 2, 3, 4, 5-point correlations. 

Although the information needed to show that the sums arising in ^ from these 
two sets of weights are the same is implicit in [5] , and although it would be easy to 
check the result on a computer, it is interesting to see what lies behind this. 
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Given a collection of disjoint subsets Ai, . . . , Am of Z/NZ and a set of weights 
Wm, k = 1, . . . m, we define c — (cq, . . . , cjv-i), where the Cj £ Z, using (jl]). The 
corresponding weighted Dirac comb is 

N-l 

and its Fourier transform is D given by 

JV-l 



S(fc)=5:c,exp('^^Up(^'=) 



where w := g-ZirV^ and P = P"^ is defined by 

The pattern frequency of (^i, . . . ,Z„) in the weighted periodic point set deter- 
mined by our choice of the cut and project scheme ([3]), the windows Ai,. . .Am, 
and the weighting system c is, up to the appropriate normalization factor, given by 



(12) Mnih, . . J„) = X CiQ+ii ■ ■ • Q+i„ ■ 

1=0 

In this way we have a function 

Af„ : (Z/iVZ)" — > R . 
A straightforward calculation of the Fourier transform of M leads to 

M^{ki,...,kn) = D{ki)...D{kn)D{-{ki+--- + kn)) 

(13) = P{w''')...P{w''")P{w-^''''+-+''-^). 

Since M„ and M„ deterimine each other, knowing one is the same as knowing 
the other. In particular, if two weighting systems determine c and c', and these 
determine the same functions P{w''^) . . . P{w^")P{w~''^^^ vkr,)^ Jq^. gQj-^g then 
they also produce the same n + 1-point pattern frequencies. 

Let us apply this analysis to the two weighting systems for Z/6Z given in pTjl . 
Here the two corresponding polynomials are 

Pi{x) := (x + l){x^ +X + l)(2x^ + 5)(3x + 1) =: Qi(x)(3a; + 1) 

and 

P2{x) := {x + l){x^+x + l){2x^ +X + 4)(3a; + 1) =: Q2{x){3x + 1) , 

Since we are only interested in the values of these at powers of w = e~^'^*/^ we may 
alter Pi and P2 so that their exponents of x are reduced modulo 6, so that they all 
lie in the range 0, . . . , 5. Having done this one checks that these indeed produce the 
coefficients cq, . . . , C5 given by wsi,ws2. We also observe that by the construction 
these polynomials vanish at w'^ ,w'^ ,w'^. Also the two polynomials are equal for x = 
1 and furthermore, w^Qi{w^^) — Q2[w) and w'^Q2{w~^') — Qi{w). This means 
that Qi{w)Piiw-^) = Q2{w)P2{w-^) and hence Pilw)Pi{w-^) = P2(w)-f2(w"^)- 
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Now consider the situation of the 5-point correlations. We have to look at 
(14) P{w^)P{w'')P{w^)P{w'-)P{w-^^+''+^+'-^) 

for the two polynomials and show that these values are the same for all possible 
values of j,k,l,r. Assuming that none of these indices is mod 6 and taking 
into account the vanishing properties above and the symmetry of the four indices 
j, k, I, r, we have the following table of possibilities: 



j 


k 


1 


r 


-(j+k+l+r) 


1 


1 


1 


1 


-4 


1 


1 


1 


-1 


-2 


1 


1 


-1 


-1 





1 


-1 


-1 


-1 


2 


-1 


-1 


-1 


-1 


4 



Given the vanishing properties of Pi, P2 , only the middle case is non-trivial. 
However, there the facts that 1 and —1 occur equally often and that the polynomials 
are equal at w = 1 gives the result. 

If one of the indices, say r is zero, then we are in the situation of the 4-point 
correlation. There the only non-obvious case is when k = ±1 and all three are 
not equal. Then along with —(i+j + k) we have 1 twice and —1 twice and so again 
the products are equal. 

The situation for the 2- and 3-point correlations is equally simple. 

We thus have: 

Theorem 2. Given the cut and project scheme ^ and a nonempty window W C 
H, then the two (different!) aperiodic weighted model sets arising from the cut and 
project scheme ([7]) and the two mod 6 weighted colourings given by wsi and WS2 
([TT|) have the same 2, 3, 4, 5-point correlations. □ 

5.1. A 2-dimensional example. The STS tiling, or shield tiling, is an aperiodic 
substitution tiling discovered by F. Gahler [5]. It consists of three types of tiles: 
squares, triangles and asymmetric hexagons looking like shields, see Fig. [TJ The 
vertices of an STS tihng can be realized as a model set .S = (R^ x , £) whose 
lattice can be described as 

C = {{ki, k2, ks, ki) : 2ki E Z for all i, all of them even or all of them odd.} . 

This lattice has a automorphisms C of order 12. In fact. Lie theorists will 
recognize the lattice as the root lattice of type F4 and C can be chosen as any one 
of its (conjugate) Coxeter transformations [2]. The four eigenvectors of C lead to 
two real C-invariant spaces and it is the projections onto these that create the cut 
and project scheme. In each of these C appears as a rotation of order 12. The 
window is a regular dodecagon W, displaced generically to avoid the projections of 
any of the lattice points of C falling on the boundary of Vl^ [8] . 

Inside C we have the sublattice of index 2 obtained by restricting the vectors 
to have integral components (the D4 lattice). Inside that there is a sublattice £0 
consisting of those vectors the sum of whose components is congruent to modulo 
3. Then C/Cq ~ Z/6Z provides us with a homomorphism a : C — > Z/6Z with 
which we can carry out the construction of fJH 

The resulting colouring on the shield tiling is indicated in the Fig. 1 with the 
different colours indicated by different symbols. Symbols that differ only by the 
presence or absence of a centre dot correspond to colours differing by 3 modulo 6. 
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Figure 1 . A fragment of the shield tihng with a 6-colouring 



There are no 6-colourings that respect the rotation C . Instead, what one can see 
here is the bands of hke-colours moving in roughly a north-by-northeast direction 
(the shortest edge vector in this direction has degree 0) and the cycling through the 
six bands as one moves in the normal directions. According to the theory above, 
each of the two weighting systems of ([TT]) can be applied to this model set and the 
resulting sets will be indistinguishable from the point of view of their 2, 3, 4, 5-point 
correlations. The results of [5] imply that they are distinguishable by their 6-point 
correlations. 

The results of on which the construction of this paper depends, seem not to 
have been generalized to dimensions greater than one. It would be interesting to 
do this since it would probably give rise to other even more interesting examples of 
distinct model sets with many identical correlations. 

6. A STOCHASTIC INTERPRETATION 

We finish by pointing out that it is possible to place the results and examples 
described in this paper into a stochastic setting. We use the same ingredients as 
before: an unweighted model set A, a homomorphism a, and a set of weights. We 
assume now that all the weights are non-negative and scaled so that they all lie in 
the range [0, 1]. We imagine now that the points of the basic unweighted model 
set A are selected or not selected on the basis of independent random choices at 
each site, the probability of being selected being Wj if the point is of colour j. The 
resulting structure is a point process, each event being the outcome of independent 



10 



XINGHUA DENG AND ROBERT V. MOODY 



Bernoulli trials made at every one of the sites according to the probabilities given 
by the weights. The moment measures that we have described then describe the 
expected values of the patterns can occur in the point process. The consequences of 
this for the two-point correlation and the corresponding diffraction measure can be 
explicitly determined from [l],Thm. 2. In particular the diffraction consists, almost 
surely, of a pure point part plus a continuous constant background. 
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